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Abstract 


In this article we prove new results regarding the existence of Bern¬ 
stein processes associated with the Cauchy problem of certain forward- 
backward systems of decoupled linear deterministic parabolic equations 
defined in Euclidean space of arbitrary dimension N £ N^, whose initial 
and final conditions are positive measures. We concentrate primarily on 
the case where the elliptic part of the parabolic operator is related to 
the Hamiltonian of an isotropic system of quantum harmonic oscillators. 
In this situation there are many Gaussian processes of interest whose 
existence follows from our analysis, including A-dimensional stationary 
and non-stationary Ornstein-Uhlenbeck processes, as well as a Bernstein 
bridge which may be interpreted as a Markovian loop in a particular 
case. We also introduce a new class of stationary non-Markovian pro¬ 
cesses which we eventually relate to the A-dimensional periodic Ornstein- 
Uhlenbeck process, and which is generated by a one-parameter family 
of non-Markovian probability measures. In this case our construction 
requires the consideration of an inhnite hierarchy of pairs of forward- 
backward heat equations associated with the pure point spectrum of the 
elliptic part, rather than just one pair in the Markovian case. We fi¬ 
nally stress the potential relevance of these new processes to statistical 
mechanics, the random evolution of loops and general pattern theory. 


1 Introduction and outline 


Let us consider the two adjoint parabolic Cauchy problems 



m(x, 0) = yoW) X G R 


N 


( 1 ) 


and 



( 2 ) 


1 


where T S (0, +oo) is arbitrary. In these equations, Ax denotes Laplace’s op¬ 
erator with respect to the spatial variable, 1^ is a real-valued function while yig 
and ijjrp are positive measures on . Both © and ([U can then be looked upon 
as defining a forward-backward system of decoupled linear deterministic heat 
equations in Euclidean space. To wit, the potential solutions to wander off to 
the future, whereas those of @ evolve into the past. In Section 2 below we show 
that under very general conditions on V,(pQ and ipj’, we can associate with (fT|)- 
m a class of the so-called Bernstein or reciprocal processes, henceforth denoted 
by These are processes that constitute a generalization of Markov 

processes which have played an increasingly important role in various areas of 
mathematics and mathematical physics over the years (see, e.g., n, m-m, cni, 
[18], [21] and the many references therein for a history and ealier works on the 
subject). In the case of ©-(El) the state space of the processes is the entire 
Euclidean space, and their construction requires a transition function as well 
as a joint probability distribution for Zq and Zt, which we denote respectively 
by P and ^ in the sequel. While P depends exclusively on the parabolic Green 
function associated with ©-©, the measure /i involves both Green’s function 
and ipQ^iprp^ or more generally a statistical mixture of (/?g and tprp. It is P and 
/i that allow one eventually to write out all the finite-dimensional distributions 
of .Z.re[o,T]) on which the remaining part of this article is based. In Section 3 
we concentrate primarily on certain Markovian Gaussian Bernstein processes 
associated with a particular case of ©-©, namely, 

1 

dtu{x, t) = - Axm(x, t)- y X (0, T], 

m(x, 0) = (/5o(x), X e (3) 

and 

1 

-dtv{x, t) = - Ax'(;(x, t)- Y 1^1^ ^)> t) e X [0, T ), 

v{x,T) = V’t(x), X e R^ (4) 

where A > 0 and |.| denotes the Euclidean norm in R^, whose right-hand side 
is the Hamiltonian of an isotropic system of quantum harmonic oscillators, up 
to a sign (see, e.g., m)- There we apply the results of Section 2 in order 
to construct and analyze three processes of interest, by determining explicitly 
in each case the finite-dimensional projections of the corresponding Gaussian 
measures. The first one is stationary and thereby esssentially a 7V-dimensional 
Ornstein-Uhlenbeck process, while the second one is a non-stationary Ornstein- 
Uhlenbeck process conditioned to start at the origin of R^. The third one is a 
process that shares many of the properties of a bridge, which we call a Bernstein 
bridge and which we may identify with a Markovian loop in a particular case. 
The measures /r we need for the construction of each one of these are intimately 
tied up with a very specific choice of initial-final data in ©-©. The situation is 
quite different in Section 4, where we construct a new family of stationary non- 
Markovian Bernstein processes that are related to the iV-dimensional periodic 
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Ornstein-Uhlenbeck process. There, we prove that the relevant non-Markovian 
probability measures /i disintegrate into statistical mixtures of the form 

Pm > 0, Pm = 1, (5) 

meN'^ meN'V 

where each is a measure related to initial-final conditions q and in 

1 

dtu{x, t) = - Axm(x, t)- y X (0, T ], 

m(x, 0) = (/j„,o(x), xeM^ (6) 

and 

1 

-dtv(x, t) = -Kv{xi, t)- Y 1^1^ ^)> t) e X [0, T), 

v{x,T) =iP^t{x), xsR^, (7) 

respectively. In other words, we show there that the construction of the ^’s 
given by ([5]) requires the consideration of a hierarchy of infinitely many pairs 
of problems of the form (O-® associated with the whole pure point spectrum 
of the elliptic operator on the right-hand side, rather than just one pair in the 
Markovian case. Finally, we also point out the potential applications of those 
new processes to statistical mechanics, to the problem of random evolution of 
loops in space and to general pattern theory, to name only three. In an appendix 
we also prove an important series expansion for the Green function associated 
with®-®. 

2 A class of Bernstein processes in 

Generally speaking Bernstein processes can take values in any topological space 
countable at infinity, and there are several equivalent ways to characterize them 
(see, e.g., [IQ]). However, the following definition will be sufficient for our pur¬ 
poses: 

Definition 1. Let N S N'*' and T S (0,-|-oo) be arbitrary. We say the 
R-'^-valued process ZT-g[o^T] defined on the complete probability space (H, J",P) 
is a Bernstein process if 

E {f{Zr) IJ-+ V JT) = E {f{Zr) \Z,,Z,) (8) 

for every bounded Borel measurable function / : R^ i—>■ R, and for all r, s, t 
satisfying r G (s,t) C [0,r]. In ®, the cr-algebras are 


7-+ = a {Z-^ {E):t<s, Eg Bn} 

(9) 

7r = CT {Z-i (7) : r > t. Eg Bn} , 

(10) 
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where Bn stands for the Borel cr-algebra on 

The dynamics of such a process at any time r G (s,t) are, therefore, solely 
determined by the properties of the process at times s and t, irrespective of its 
behavior prior to instant s and after instant t. Of course, it is this fact that 
generalizes the Markov property. 

In order to associate a mere Bernstein process with m-m we now impose 
the following hypothesis, which regards Green’s function alone: 

(HI) The measurable function V : R is such that the parabolic 

Green function g associated with ©-© is jointly continuous in all variables 
and satisfies 


5 (x,t,y) > 0 

for all X, y G R^ and every t G (0, T], 


( 11 ) 


Let us now write Ai (R^ x R^,C) for the space of measures we are inter¬ 
ested in, namely, the topological dual of the Frechet space of all complex-valued, 
compactly supported functions on R^ x R^ endowed with the usual locally con¬ 
vex topology (see, e.g., [H]). Having (fTTl) at our disposal, let us introduce the 
functions 


p(x,t;z,r;y,s) 


g(x,t-r, z)g(z,r-s,y) 
g{x,t- s,y) 


( 12 ) 


and 


P{x,t;E,r-y,s) 


j dzp (x, t; z, r; y, s) 

E 


(13) 


for every E G Bn, both being well defined and positive for all x,y, z G R'^ and 
all r,s,t satisfying r G {s,t) C [0,T]. For every F G Bn x Bn, let us also 
consider a positive measure g, G M (M'^ x R^, C) such that 


g{F) := [ dg{x,y) (14) 

JF 

defines a probability measure on x Bn, thus satisfying 

f dg{x,y) = l. (15) 

The knowledge of both (US and (HI then makes it possible to associate a 
Bernstein process with (HD-©. The precise statement is the following: 


Theorem 1. Assume that Hypothesis (HI) holds, and let P and g be 
given by m and (B-m, respectively. Then, there exists a probability space 
(H,J^, P^) supporting an -valued Bernstein process .Z,-g[o^T] such that the 
following properties are valid: 
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(a) The function P is the transition function of in the sense that 

{Zr G E \Zg, Zt) = P {Z-t, t] E, r; Zg, s) 

for each E € Bn and all r,s,t satisfying r € {s,t) C [0,T]. 

(b) For every n S N'*', the finite-dimensional distributions of the process are 
given by 


IP/i (-^ti G El,Zt„ € En) 

r dn{x,y) r 
Jr^xs.^ 5(x, T’, y) Jei 


dx„ 


n 

X 5(xfe,ifc - 4-i,Xfe_i) X 5(y,r-t„,x„) 

k=l 


(16) 


for all El, .■.,En € Bn and all to = 0 < ti < ... < tn < T, where Xg = x. In 
particular we have 


P/. [Zt e E) 



f dfi[x,y) 

[ dzg[K,t,z)g[z,T -t,y) 

(17) 


JrNxR'v 5(,x, a ,y) , 

J E 

for each E G Bn and every t G (0,T). 

Moreover, 



P/i [Zo G E) 

= Pl[Ex R^) 

(18) 

and 





Pm [Zt e E) 

X 

A- 

a. 

II 

(19) 


for each E G Bn. 

(c) IPis the only probability measure leading to the above properties. 


Proof. Up to minor technical details, a direct adaptation of the method 
developed in Section 2 of [21] leads to 


P/i [Zo & Eq, Zt„ G El,..., Zt.^ G En, Zt G Et) 



d'Xn W_p[y,T-,Xk,tk-,Xk-i,tk-i) 

k=l 


( 20 ) 


for all Eo, ..■,Et G Bn and all tg = 0 < < ... < < T, where Xg = x. In 

particular we have 


P/i [Zo G Eo, Zt G Et) = /r [Eo x Et) (21) 

for all Eo, Et & Bn, that is, m is the joint probability distribution of Zo and 
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Zt. Now, from (HU we obtain 


n 

P (Vi i^k-) ^k—l ) —1) 


n 


k^l 


5f(y,T-4_i,Xfc-i) 


1 

9i>^,T,y) 


n 

JJ p(xfc,4 - 4-i,Xfc_i) X p(y,T - tn,x^) 

fe=i 


after n — 1 cancellations of factors in order to obtain the second equality, so that 
(fTCl) follows by choosing Eq = Et = in (pHl) . Relations (fTSl) and (|TO1) are 
just a particular case of (EU. ■ 


Of course, we can say more about if we know more about /i. First, 

^Te[o,T] is Markovian if there exist positive measures vq, vt & M. (R.^, C) such 
that 

fi{F)=[ O i/t) (x,y)g(x,T,y) (22) 

J F 

for every F G Bn Bn, with 

[ d(z/o ® i^t) (x,y)5(x,'7",y) = 1 (23) 

in which case we also say that p, is Markovian. This result can be traced back to 
the more general Theorem 3.1 in uni, and allows us to make a closer connection 
between Zt-^[o^t] and (H])-© provided we impose the following hypothesis: 

(H2) The measures (Pq,iPt ^ (M^,C) in 0-0 are positive, and there 

exist a unique classical positive solution to o and a unique classical positive 
solution to E]), namely, 


Uvoi^,t)= d(po{y)g{x,t,y) (24) 

and 

v^^{x,t)=[ dN{y) 9 ix,T-t,y), (25) 

JRW 

respectively. 

We then have the following consequence of Theorem 1: 

Corollary 1. Assume that Hypotheses (H1)-(H2) hold, and let us choose 
pL of the form I122\} with vq = and vt = '0r- Then, the Bernstein process 
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Zt^[o,t] of Theorem 1 is Markovian, and its finite-dimensional distributions are 
given by 


P/i (-^ti G El,Zt„ € En) (26) 

I dxi... I dxji I I g (xfc ,t[„ tf^—i, ^ u^p^ (x^, ti'jv.ip.^ (xj^, tn) 

J JE„ o 


k=2 


for every integer n> 2, all Ei,..., En € Bn and all tg = 0 < ti < ... < tn < T, 
where Xp = x. Moreover we have 


Pm {Zt G E) 

/ dxup^{x,t)v^.^{x,t) (27) 


for each E G Bn and every t G (0,T). Finally, 



Pm {Zo & E) — 

/ d>pQ (x) (x, 0) 

J E 

(28) 

and 

P^ {Zt GE) = 

/ (x) (x, T) 

J E 

(29) 

for each E & Bn- 





Proof. We first rewrite (fTOl) as 


P/i {Zti G El,..., Zt 
f dfi{x,y) r 

/r" xR-^ 5(x,'P,y) Jei 


G En 
dxi. 


dXn 


I En 


n 

X n 5(xfe,tfe - 4-i,Xfc_i) X g {xi,ti,x) g {y,T - tn,Xn) 
k=2 


(30) 


and then substitute the choice of the measure /i into dUD and (I30L using (1241) - 
(|25l) along with the symmetry property of g with respect to the spatial variables. 


Remarks. (1) In the Markovian case, we have thus exhibited a very gen¬ 
eral class of initial-final conditions in 0-0 which allows us to determine the 
processes Z.re[o.T] completely, including their marginal distributions (l^ - (l^ . 
The converse point of view was developed in [5] and its references, where it was 
shown instead that it is a general class of marginal distributions which deter¬ 
mines the initial-final data of the relevant partial differential equations, through 
a system of nonlinear integral equations. However, the resulting initial-final 
conditions of [2] belong to the class of positive continuous functions, and not to 
the larger class of positive measures as is the case in this article. 
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(2) A glance at (1^51) shows that it is sufficient to use elementary time reversal 
in Green’s function to obtain the solution to ([2]). Although the situation is not 
always that simple, particularly when the given parabolic equations are non- 
autonomous, it is still possible to define a quite appropriate probabilistic notion 
of time symmetry in general. We refer the reader to |21] for further details. 

Relations dm) and (l26t are the fundamental relations that will allow us to 
construct the Gaussian processes associated with ©-(HI) in the next sections. 


3 Two Ornstein-Uhlenbeck processes and a Bern¬ 
stein bridge 


In the remaining part of this article we denote by (., .)]jjv the Euclidean inner 
product in , and by the usual Lebesgue space of all complex¬ 

valued, square-integrable functions on . We begin by considering the forward- 
backward system ©-([I]) with centered Gaussian initial-final data, namely. 


dtu{x, t) 
u(x, 0) 


iAxu(x,t) 
‘/5o,aW = 



(x,t) G R^x (0,T], 


A exp [AT] 


exp 


A|xr 

2 


(31) 


and 

- dtv(x, t) 

u(x,r) 


^A^v(x,t) - ^ u(x,t), (x,t) G R^x [0,r), 




/ A exp [AT] 


exp 


A lx 


(32) 


thereby identifying the measure a = V’t a with its Gaussian density relative 
to the Lebesgue measure in R^. Let us recall that the self-adjoint realization 
in (R^,C) of the elliptic operator on the right-hand side of these equations 
has a pure point spectrum. More specifically, for every m G N let 

hm,\{x) := vT/im (^a/Ax^ (33) 

be the one-dimensional, suitably scaled Hermite functions where 

hm {x) = (^7r52'"TO!^ e~^Hm{x), (34) 

and where the iLm’s are the Hermite polynomials 

2 dr 

dx^ 


Hm{x) = (-l)™e"' 


(35) 













The following spectral result regarding the operator on the right-hand side of 
ED is well-known and in any case can be verified directly by means of an explicit 
computation; 


Lemma 1. The tensor products where the rrij’s vary indepen¬ 

dently on N provide an orthonormal basis of (]R^,C), and moreover the 
eigenvalue equation 



f=l hmj,\ (x) = A 



+ T I ®f=i 


,.a(x) (36) 


holds for every x gK.^. 


In particular, by reference to (IMl) and (1551) we have 


li /lo.A (x) = 


N 

1- 

CM 

1_ 

- exp 

— 

\TrJ 

2 


> 0 


(37) 


whose associated eigenvalue is E\ = ^, so that (l37l) corresponds to the initial- 
final conditions in (I31l) - (l32p up to a normalization factor chosen in such a way 
that m holds for the measure 


Tx (F) = / dxdyipo^), (x) V't.a (y) 5a(x, T, y) 

JF 

where F G Bn x Bn, and where 

gx{x,t,y) (38) 


( ^ ) 

K 

exp 

A 

cosh(At) 

(|x|V|y|^) 

2 (x, y)]Riv^ 

\^27r sinh (At) ^ 



2 sinh (At) 



is the iV-dimensional version of Mehler’s kernel for t G (0,T] (see the appendix 
for details). In this way, the unique classical positive solutions to (I^ - (I5^ 
satisfying the requirement 


/ dxuA(x,f)uA(x,t) = 1 


for every t G [0, T] are 


, , /Aexp [ATI 

UA(x,t) = ( ---^ ) exp 


A ixr +ivt 


and 


, , / A exp [AT] 

vxix,t) = ( - - -1 ) exp 


Afjxr+ fV(T-t 


(39) 


(40) 


(41) 
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respectively. Then the following result holds, where stands for the expecta¬ 
tion functional on (17, P^): 


Proposition 1. The Bernstein process 7 ^] associated with m-m 

in the sense of Corollary 1 is a Gaussian and Markovian process such that 


£ E) = {2TTa\) ^ / dxexp 
J E 


|x| 


for each t € [0,T] and every E € Bn, where 

1 


cta = 


2A' 


Furthermore, the components of satisfy the relation 




exp [—A \t — s|] 

— 2 ^ 




(42) 


(43) 


for all s,t G [0, T] and all i,j G {1, N}. Thus, Z^^^^ identifies in law with 

a process whose components are all independent, one-dimensional and stationary 
Ornstein- Uhlenbeck processes. 


Proof. The fact that Z^^^^ is Gaussian and Markovian satisfying (I42p 
follows from Corollary 1 with (1^ and (|i(I)) - (HT]) plugged into (1^ and (l?7l) . so 
that the density of the probability distribution for [Z^^, ..., Zf^^ G is 


9\ {^k,tk - tfc-l,Xfc_i) X UA(xi,ti)uA(x„,t„) 


/ n > 

= ( 2 ^)-^ ( 2 A"e^(‘"-‘i)) ' n 


\k^2 


X exp 


^ cosh(A(tfe - tfe_i)) (|xfc| -f|Xfc_i| j- 

~2 ^ 
fc =2 


sinh (A(tfe - tfc-i)) 


K 

2 


2 {f■k^'>^k-^) 


R" 


X exp 



The nN x nN corresponding covariance matrix is then of the form Ca 0 Hat 
with In the identity matrix in K^, where is tridiagonal and obtained by 
identification of the quadratic form in the argument of the above exponentials 
when TV = 1. This gives 


Cxik = 


Xe^G2-tl) 

sinh(A(t2—ti)) 
A sinh(A(t|c+i-t|i;i)) 


for k = I, 


sinh(A(ifc+i-tfc)) sinh(A(tfc-tfc_i)) 


for k = 2, ..., n — 1 , 


sinh(A(tn —in -1)) 


for k = n 
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(the second line not being there when n = 2), and 


1 _ 


_A_ 

sinh(A |4 - 4-11) 


for k = 2,n. 


Consequently, using for instance the analytical inversion formulae in [9], or by 
direct verification, we obtain 


exp[-A|4-til] 

- 

for all k,l G {1,..., n}, so that (l4^ eventually holds. Now, let us consider the 
forward Ito integral equation 




^0 


„ — Xt r^X I 

e Zp + 
^0 


Jo 


tG [0,T], 


(44) 


where Wr^[o,T] is a given Wiener process in R^, and where Zq is distributed 
according to (142 p and independent of ■ It is well known that the solution 

process r] to (l44l) is centered Gaussian with covariance (|^ (see, e.g.. 

Section 5.6 in Chapter 5 of [II]), so that Z^^jp identifies in law with ^^g[p yj. 
Therefore, Z^^ jp yj is indeed a A^-dimensional Ornstein-Uhlenbeck process with 
the stated properties. ■ 


Next, we show that if we require instead Zp to be a given point in K'^, the 
solution process to (I44p is also a particular Bernstein process. We do this in 
the simplest case where Zp is the origin, by considering the forward-backward 
system 


and 


1 A^ Ixl^ 

dtu{x,t) = -Axu(x,<)-^^(x,^), (x,t) G K^x (0,r], 

m(x, 0) = (/3p,a(x) = (exp [AT])^ 5(x) 


(45) 


-dtv{x,t) = ^A^v{x,t) - ^ v{x,t), 

v(x,T) = = (exp [AT])^ exp 


(x,t) G R^x [o,r), 

Alx^ ' 

2 


(46) 


with S the Dirac measure. The corresponding Markovian measure (HI is then 
determined by 


AiA(x.y) = (exp [AT]) 2 <5(x)exp 


A|y|^ 

2 


5(x, T,y) 
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and satisfies m, while the unique relevant classical positive solutions to (HSl) - 
((46l) are 


ux{x,t) 


fwxirexp 

^27r sinh(At) j 


Acoth(At) |x|^ 
2 


(47) 


and 


vx{>^,t) = (exp [AT]) * exp 


A(|x|2 + iV(T-t)) 
2 


(48) 


respectively. We then have the following result: 


Proposition 2. The Bernstein process associated with |^5[ )-|.^6|j 

in the sense of Corollary 1 is a Gaussian and Markovian process such that 


P;, G E) = (27rcrx(t)) 2 f dxexp 

J E 

for each t G (0,T] and every E G Bn, where 

sinh(At) exp [—At] 

T-• 


2crA(t) 


Furthermore we have 


(4 = o) = 1, 


and the components of satisfy the relation 


E,; 


+ (exp [2A it A s)] - 1) <5,,,- 


(49) 


(50) 

(51) 


(52) 


for all s,t G [0,T] and all i,j G {1, Thus, Z^^^^ identifies in law with a 

process whose components are all independent, one-dimensional non-stationary 
Ornstein- Uhlenbeck processes conditioned to start at the origin. 


Proof. As in the proof of Proposition 1, the first part of the statement 
including (l4^ - (ICTl follows from the appropriate substitutions into the formulae 
of Corollary 1. Furthermore, the matrix resulting from the identification 
of the quadratic form in the Gaussian density of iZf'^,Zf'^) is the same as 
that in Proposition 1, with the exception of which now reads 

__ Asinh(At2) _ 

sinh (A(t 2 — ti)) sinh (Ati) 

Inverting again we eventually obtain 

c**. = („p |2A (,, A i,)l -1) 
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for all k,l S n}, so that (15^ holds. Therefore, is indeed a N- 

dimensional Ornstein-Uhlenbeck process with the stated properties. ■ 


Remarks. (1) It is equally easy to condition the Ornstein-Uhlenbeck process 
so that it ends at the origin of when t = T. The underlying Bernstein process 
^Te[o T] then simply determined by swaping the initial-final conditions in (1451) 
and (gSl). Indeed, in doing so the relevant solutions become 


K 

u\{K,t) = (exp [AT]) ^ exp 


A 



2 


and 

Acoth(A(T — t)) |x|^ 

2 

instead of (l47l) and (l48l) . respectively, so that is just the time-reversal 

of the process of Proposition 2, namely, 

for every re [0, T]. 

(2) Whereas the Bernstein process of Proposition 1 is stationary, that of 
Proposition 2 is not. This is intuitively understandable, as some kind of non 
trivial dynamics ought to be necessary to steer the process from a deterministic 
state at t = 0 to a Gaussian distribution at the end of its journey. 


vx{x,t) = 


A exp [ 




27rsinh(A(r — t)) 


exp 


Of course, there is an almost unlimited number of possibilities of getting 
various Bernstein processes in the above manner, just by choosing cpQ ^ and ijjrp ^ 
appropriately. We complete this section by constructing yet another process 
which shares many properties of a Markovian bridge. For this, we consider 
the forward-backward system dSD-O with Dirac measures concentrated at the 
origin and at a given point a G as initial-final data, namely, 


A" IxP 


dtu{x, t) = - Axu(x, t) -t) e R X (0, T], 

u(x,0) = (^o.aM = mA5(x) 


(53) 


and 


A' IxP 


-dtv{x,t) =-A^v(x,t) - - — v(x,t), (x,t)GR x[0,r), 

u(x, T) = '0 t,a(x) = mA<5(x - a), 


( 54 ) 


where 


rriA = 


27r sinh 




^ / 


exp 


Acoth(AT) |a| 
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It is easily verified that the normalization condition 


m. 


xR^ 


dxdy6{x)S{y - a)gx{x,T,y) = 1 


holds with Mehler’s kernel given above, and that the unique classical positive 
solutions to (l53l) - (l54l) which satisfy (l3^ are now 


itA(x, i) = riA sinh 2 (At) exp 


axi,t) |x|' 


(55) 


and 


vx{x,t) = Ha exp 


axiT - t) |a| 


sinh-^ {\{T-t)) 


X exp 


-- aA(T-t)|x|^- 


2A (a,x) 


sinh (A (T — t)) J 


respectively, where 
and 


ax{t) = Acoth(At) 


Ha = 


N 

/Asinh(AT)\T 

UA(T)|a|^l 

(, 2;r ) ""P 

4 


(56) 

(57) 

(58) 

(59) 


Then the following result is valid: 

Proposition 3. The Bernstein process associated with {50- 154]) in 

the sense of Corollary 1 is a Gaussian and Markovian process such that 


{z^ G E) = {2TTax{t)) ” [ dxexp 

J E 

for each t G (0,T) and every E G Bn, where 


|x- aA(t)r 

2ax{t) 


and 


Furthermore we have 




, , sinh(At) 

sinh (A(T — t)) sinh(At) 


A sinh (AT) 


\ (^0 = o ) 

(^T = a)=l, 


(60) 

(61) 

(62) 


(63) 
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and the components of x] satisfy the relation 


E, 


- a\is)){Z^-^ - a 


sinh(A(T—t)) sinh(As) 
Asinh(AT) ' 




for t>s, 


sinh(A(T-s))sinh(At) c fnr ^ « 

AiSh(AT) tor t S S 


(64) 

for all s,t G [0, T] and all i,j € {1, N}. In fact, Z^^ jQ yj is a non-stationary 
process pinned down at the origin when t = 0, at 3 when t = T, and exhibiting 
maximal randomness when t = ^. 

Proof. We begin by proving (1501) . Using (1551) . (1551) and (1551) we first have 

ux{x,t)vx{x,t) 


X sinh(AT) 


-V 


27rsinh (A(T — t)) sinh(At) J 
1 


exp 


{ax{T)-ax{T-t))\B\ 


“2 ( 1^1 “77 


X exp 

after regrouping terms, and furthermore 

a\{T) - ax{T - t) = - — 

axit)axiT - t) = 


2A (a, x)][jAr 


sinh (A (T — t)) 

A sinh(At) 
sinh (A(T — t)) sinh(Ar) 
A sinh (AT) 


sinh (A(r — t)) sinh(At) 
from ((58l) . The substitution of (l66l) - (l67)) into (l65)) then leads to 
ux{^,t)vx{xi,t) 


(65) 

( 66 ) 

(67) 


A sinh(AT) 


exp 


Asinh(At) |a| 


27rsinh (A(T — t)) sinh(At) 

A /sinh(AT) |x|^ — 2sinh(At) (a,x 
2 y sinh (A(T — t)) sinh(At) 


2 sinh (A(T — t)) sinh(AT) 


X exp 


Now, for the numerator of the argument in the second exponential of 
have 


( 68 ) 


we 


smh(AT) \x\ — 2sinh(At) {a,x)^N 

. , , X m 2 sinh^(At) |a| 

= sinh(AT)|x-aA(t)| - ^ ' 


sinh (AT) 

by virtue of ((5T1) . Therefore, taking (1551) and into account in 

uxi>^,t)vx{^,t) 


(69) 
we get 


/ XX — K 

= {2TTax{t)) 2 exp 


|x- aA(0l^ 

2cjx{t) 
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following the cancellation of two exponential factors, which gives the desired 
result according to (1271) . 

As for the proof of (I55)) . we remark that (175)) and (1771) lead to 


P/j {Zo S E) 


f 

A ^cosh(Ar) X ^ — 2 (a,x)][jN^ 

1 LL/\(J 1 

J E 

2 sinh (AT) 


and 


{Zt £ E) 


[«A(T)|a|2l 


\ «A(T)|xn 

2 

/ dxd[x — a)exp 

Je 

2 


for every E £ Bn, respectively, which immediately imply the claim. 

We now turn to the proof of (1641) . by noticing that in this case the finite¬ 
dimensional density in is 


gx {^k,tk - tfc-l,Xfe_i) X MA(xi,tl)uA(x„,t„) 

fc =2 


= (27r) ^ 


A” sinh(AT) 


sinh (A(T — t„)) sinh(Ati) 


sinh(A(4 - tk-i)) 


\k^2 


X exp 


i(aA(T) - ax{T-tn)) |a|^ 


H 

2 


X exp 

X exp 


^ cosh(A(tfc - 4_i)) (|Xfc|^ -k |Xfe_i|^) - 2 (xfc,Xfc_i)jg,^ 
“o' 


fe =2 


sinh(A(tfe - tk-i)) 


^ (a\{ti) |xi|^ -I- ax(T - t„) |x„|^^ 


X exp 


A (a,x„ 


sinh(A(r - tn)) 


Therefore, for the tridiagonal matrix ^ corresponding to the quadratic part 
when = 1 we obtain 


Asinh(At 2 ) r. — 1 

sinh(A(t 2 —ii)) sinh(Ati) ’ 


Cxik = 


_Asinh(A(tfc+i-t|i;_i))_ 

sinh(A(tj 5 +i-tfc)) sinh(A(tfe-tfc_i)) 

Asinh(A(T-t„i)) 


sinh(A(T—in)) sinh(A(in—in-l)) 


for fc = 2,..., n — 1, 
for k = n 


(the second line still not being there if n = 2), and 


1 _ 


_A_ 

sinh (A \tk-tk-i\) 


for k = 2 ,..., n. 
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Consequently, inverting again and using numerous relations among hyperbolic 
functions we eventually get 

{ sinh(A(T-tfc))sinh(Ati) fny. > / 

Asinh(AT) lOr K ^ t, 

sinh(A(r-ti))sinh(Atfc) I, < ; 

sinh( AT) K S b 

which leads to (IMl) by standard arguments. Finally, we note that the curve a\ : 
[0,T] H->- RJ given by (1621) is concave aside from satisfying ctaIO) = <y\{T) = 0, 
and that it takes on the maximal value 

/r\ sinh^(^) 

^\2) Asinh(AT)’ 

Thus, the process yj is indeed non-stationary with the stated properties. 

■ 

Remark. We may dub the process t] preceding corollary a 

Bernstein bridge, which represents a random curve whose ends are pinned down 
at specified points in space. We remark that the corresponding Gaussian law 
is no longer centered, unless a = o in which case the process materializes a 
Markovian loop which retains the main features of a Brownian loop. In fact, 
y] does reduce to a Brownian bridge in the limit A —?> 0+ since 

5ij for t> s, 
j for t < s 

according to (l64l) . 

In the next section we introduce a new class of Bernstein processes which 
we can eventually relate to the so-called periodic Ornstein-Uhlenbeck process, 
and which is generated by a one-parameter family of non-Markovian probability 
measures. 


^hm - ai(t))) 


4 A family of non-Markovian Bernstein processes 
and the periodic Ornstein-Uhlenbeck process 

The Bernstein processes of this section are still defined from measures which are 
intimately tied up with problems of the form ®-0: but their finite-dimensional 
distributions will be determined exclusively from Theorem 1. Instead of consid¬ 
ering just one pair of parabolic problems such as ©-0, we first introduce an 
infinite hierarchy of forward-backward systems of the form 

1 A^ 

dtu{K,t) = -Axu(x,t) - y |x|^M(x,t), (x,t) e X (0,T], 

n(x,O) = (p,,_o^;,(x)=eMs:f=i-i+fK 0 f^^h™^_;,(x), x e (70) 
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and 


1 

-dtv{x, t) = - Ax-!;(x, i) - y |x|^ v[x, t), (x, <) e x [0, T), 

■u(x,r) = ^/>^ T_;^(x) = eMSf=i'"^ + f (g)|^Ym,.A(x), xeR^, (71) 

namely, one such pair for each rrij G N and every j, where the hmj,\s are the 
Hermite functions of Lemma 1. Accordingly, this means that we are considering 
as many pairs of such systems as is necessary to take into account the whole 
pure point spectrum of the elliptic operator on the right-hand side. We remark 
that GOD-dni) constitutes a generalization of (ED)-®, the latter system being 
associated with the bottom of the spectrum where rnj= 0 for each j. Whereas 
the associated measures remain suitably normalized, the drawback is that they 
are no longer positive according to the following result: 

Lemma 2. Let us write m := (mi,...,mAr) G and let us consider the 
sequence of measures 

= f C«Mm,A(x,y) (72) 

JF 

where 

Mm.A (x, y) = V^m.O.A (x) V'rn.T.A (y) 5'a(x, T, y) (73) 

and F G Bn x Bn- Then we have 

[ (x, y) = 1 (74) 

and for every m 0 the ^^.a’s signed measures. 


Proof. According to Proposition A.l of the appendix we have 
5a(x, r,y) 

= (x) X (y) (75) 


neN'^ 


for Mehler’s kernel (I38p . where the series converges absolutely and uniformly 
for all x,y G Moreover, since the tensor products a provide an 

orthonormal basis of (R^,C) we have 


N 

/ dx (x) X 0^1 A (x) = TT Smj, 


(76) 


for all rnj,nj G N. Consequently, substituting g a, '0m t a ED into the 
left-hand side of (1711) and taking (1751) into account we obtain 


/R^xR^ 


(x) 0,^ T_A (y) 5 a(x, T, y) 


N 


ngN" i=l 
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Finally, the are signed measures because of the existence of real zeroes 

for each hmj.x when mj ^0. I 


The fact that the above measures are signed prevents one from applying 
directly the general results of Section 2 to (TTOt - dTI]) with a fixed m ^ 0, as it 
would prevent one from applying the main result of [2] briefly discussed at the 
end of Section 2. However, we can still save the day by constructing a one- 
parameter family of bona fide probability measures from all the _),’s. Indeed 
the following result is valid, where we call a measure non-Markovian whenever 
((^ does not hold: 


Lemma 3. For each A > 0 there exists a one-parameter family 
of positive, non-Markovian measures of the form 

t^\AF)= / 

JF 

where F € Bn x Bn, which satisfy 




(x, y) = 1 


(77) 


and which disintegrate into a statistical mixture of the ^ dn other words, 
for each A > 0, each 6 > 0 and every m G there exist numbers Pm,x,e > 0 
such that 

d\,e = X! ^ p^.A.s = 1- 

In fact, it is sufficient to take 


AA,e (X) y) = (2 (cosh(A {6 -\-l)T)- 1)) ^ gffix, T, y)gx{x, 6»T, y). (78) 


Proof. From the series expansion (1751) and the orthogonality relations (1761) 
we have 


/ dxdygx (x, T, y)gx (x, OF, y) 
'r'^xR'v 

leN'v 

-N 


= I 2 sinh ( A (0 -I- 1) 


T 


by summing the underlying geometric series, from which we obtain dZZl) as a 
consequence of ([75|l and the identity 

dsinh^ ( A(6(-f 1)^ ) = 2 (cosh(A (d + 1) T) - 1). 
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As for the second statement of the lemma, let us define the sequence of numbers 


Pm.A.e := (2 (cosh(A (B + 1) T) - 1))^ + f )^(»+i)t > q. (79) 

Summing as above we get 

^ ^ Pm,X,& — f 
meN^ 

as required. Therefore, taking into account (Ea), (na and dzi we obtain 
Pm,A,e/^m,A (X>y) 

= 5A(x,T,y) Y Pm.A.ey^m.O.A WV’m,T,A(y) 

= (2(cosh(A(0 + l)r)-l))^<7A(x,r,y) 

X Y m, + f)\9T ^ (8>Y^rn,.X (y) 

meN« 

= AA,e(x,y) 

according to (1751) . which is the desired conclusion. ■ 


Remark. Strictly speaking, the measure does not exist for 0 = 0 but 
the limit 

Aa.+ (x,y) := lim /iA,e (x, y) = (2 (cosh(Ar) - l))"^ gA(x, T, y)(5(x - y) (80) 

C7^U-f 

does, by virtue of the fact that gx is Green’s function associated with the partial 
differential equation in (1311) . Said differently. Lemma 3 and its proof remain valid 
for the measure jix.+ associated with (1501) since we have 

Y (x) X ®Yhmj,x {y) = S{x - y) 

in the sense of distributions for all x, y as a consequence of the complete¬ 

ness of the (8)jLi^mj,A’s in (]R^,C). 


The above developments now lead to the following result: 


Proposition 4. For each A > 0, the Bernstein processes associ¬ 

ated with the above infinite hierarchy through the measures fix e stationary 
Gaussian and non-Markovian processes such that 


P, 


G = (27r(TA,e) J dxexp 


2a 


X,0 


(81) 
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for each t € [0,r] and every E G Bn, where 


crx,0 = 


sinh (A {0 + 1) T) 
2A(cosh(A(6l + l)r)- 1)' 
■\,e 


Furthermore, the components of satisfy the relation 

cosh(^x(\t-s\ - 


E, 




2Asinh(^(^) 




(82) 


(83) 


for all s,t G [0, T] and all i,j G {1, N}. Finally, the process t] ‘associ¬ 
ated with the measure _|_ identifies in law with a stationary process whose com¬ 
ponents are all independent, one-dimensional and periodic Ornstein-Uhlenbeck 
processes. 


Proof. The processes are Gaussian and non-Markovian by virtue 

of (fTBl) with Green’s function (155)) and the measures Furthermore, using 
the symmetry properties of g\ with respect to the spatial variables, and twice 
the related semi-group composition law (IMl) of the appendix, (HU with dUI) 
becomes 

G f;) 

= (2 (cosh(A (0-I-1) T) - 1))~ / dxg\{x,{9 + l)T,x) 

Je 

for every t G (0,T). The same result obtains for t = 0 and t = T as a conse¬ 
quence of (fT5)) - (fT51) . so that ((5T]1 - (I5511 follows immediately since 


5a(x, (0 + l)r,x) 

_A_ 

27r sinh (A (0 -(- 1) T) 


exp 


A (cosh(A (6 >-fl) T) - 1) |x|' 
sinh (A {6 + 1) T) 


We now turn to the proof of (j83|) by determining the Gaussian density of 

( X 0 X 0\ 

Z(^’ ,Z(^’ 1. Thus, using px^g in ()16p and integrating first over y and then 
over X we may write 

p^(zy GF;i,...,z*y g£;„) 

= (2 (cosh(A (0-I-1) T) — 1))~ f dxi... ( dx„ 




n p 

Y\ g\{xk,tk-tk-i,Xk-i) X / dxgx{xi,ti,x) gxix,{9 -\-l)T -tn,Xn) 

k=2 


= (2 (cosh(A (0-I-1) T) — 1)) 2 / dxi... / dx„ 

J Ex 'J Eji 


X n (Xfe,tfe - 4-i,Xfc_i) X gx (xi, {9-\-l)T- {tn - ti),x„) 
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where we have once again used the semi-group composition law (IMl) twice, so 
that the Gaussian density in reads 


= (2tt 


(2(cosh(A(6i-Hl)T)-l))“ 

n 

X n (Xfc,4 - 4-i,Xfc_i) X gx (xi, {9 + l)T- {tn - ti),x„) 

k=2 

f 2A"(cosh(A(6»-hl)T)-l) 


X exp 


X exp 


sinh(A((0+l)T-(t„-ti)))j n(sinMA(ife ifc-i)) 
^ cosh(A(4 - tk-i)) (|xfc|^ -h |xfc_i|^) - 2 (xfc,Xfe_i)jj„ 


k^2 


sinh(A(tfe - tk-i)) 


^ cosh{X{{0+ 1)T - {tn-h))) [\xi\ -F|x„| j - 2(xi,x„)jj«^ 
2 sinh(A((6»-I- l)r - - ti))) 


For the sake of clarity we now identify the corresponding matrix I by con¬ 
sidering the case n = 2 separately from the case n > 3, as is no longer 
tridagonal. For n = 2 we obtain 


„_i _ Asinh(A (0-I-1)T) 

x,e,k,k {\{t 2 — ti)) sinh (A((0 + 1) T — (t 2 — G))) 


for fc = 1, 2 


and 


^A,0,2.1 ~ ^A,0,1.2 


while for n > 3 we get 


A 


A 


sinh (A 1^2 - G|) sinh(A((0 -I- 1) T - 1^2 - G|)) ’ 


Asinh(A((e-|-l)T-(tn-t 2 ))) 


for k = 1, 


^\},k,k — 


sinh(A(t 2 — )) sinh(A((0+l)T— (tn — ti))) 

Asinh(Afa+,-t._i)) - forfc = 2,...,n-l, 


sinh(A(t|e+i-tfc)) sinh(A(tfc-t|e_i)) 

_ Asinh(A((e+l)T—(tn_i —ti))) _ r , 

sinh(A(t„-t„_i))sinh(A((e-|-l)T-(i„-ti))) lUl H — n , 


1 ___ 

^A,0,fe,fe-1 ~ ^A,0,/c-l,fc ~ 


-1 


and 


sinh(A |4 - 4-i|) ^ 2, 


A 


C'A.e.n.i C'A.e.i.n sinh (A((0-h 1) r - |t„-til))’ 


all the remaining matrix elements being zero. In both cases we then obtain by 
inversion 


Cx, 


0,k,l = 


sinh (A |tfe - til) - sinh (A(|tfe - ti| - (0 -I- 1) T)) 
2A (cosh(A (0 + 1) T) - 1) 
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for all k,l € {1,n} or, equivalently, 




cosh (^X (^\tk -ti\- 

2Asinh(^^^^) 


so that (l83|) eventually holds. Let us now consider the case 0 = 0, namely, the 
case corresponding to the measure flx.+ defined from ([80]), for which (1^ and 
(|Mjl become 


sinh (AT) 

2A (cosh(Ar) - 1) 


(84) 


and 


E, 




cosh (A(|t- s| - f)) 
2 Asinh(^) 


(85) 


for alH, j S {1, ■■■, N}, respectively. Let us also consider the forward Ito integral 
equation with periodic boundary conditions 




^0 


e-^*Xo+ [ t€[0,T], 

Jo 


Xx 


( 86 ) 


rather than (|44p . It is known from Theorem 2.1 in [14] (see also m or Section 5 
in [T5] for the case = 1) that the solution to (15^ can be written out explicitly 
and defines a non-Markovian centered Gaussian stationary process, namely, the 
so-called periodic Ornstein-Uhlenbeck process given by 

I ^ t e [0,T], 

whose variance and covariance are given by (IMl) and (1551) . respectively. There¬ 
fore, identifies in law with that process. I 


Remarks. (1) It is of course also a posteriori clear that the processes of 
Proposition 4 are not Markovian since the covariances (1551) do not factorize as 
the product of a function of s times a function of t, in contrast to all the cases 
investigated in Section 3. Furthermore, a very different way of understanding 
such factorization properties in the Markovian case was put forward in Section 
6 of [18] . where the covariances were written as the product of two linearly 
independent solutions to some suitable Sturm-Liouville problems. 

(2) Problem (1551) is part of a more general class of linear stochastic differential 
equations that were investigated by several authors, including [14], m-m 
and some of the references therein. In this context we ought to mention an 
analysis of the law of the solution to (1861) in one dimension carried out in Section 
5 of [15], which establishes a relation with Bernstein processes whose state 
space is one-dimensional. The main tool used there is a formula of integration 
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by parts proved directly on the underlying infinite-dimensional path space by 
means of Malliavin’s calculus. This is in sharp contrast to the method we have 
developed in this section, as we have first constructed a one-parameter family of 
non-Markovian, A^-dimensional Bernstein processes associated with the infinite 
system GnU-dllD, which we have only a posteriori identified with the solution 
to (155)) when 9 = 0. 

(3) While the processes materialize a one-parameter family of non- 

Markovian random curves in , they might also occur naturally in completely 
different contexts, as the limiting process Zj^^^ does. Thus, this process is 
quite relevant to the mathematical investigation of certain quantum systems in 
equilibrium with a thermal bath since it identifies with the Gaussian process 
of mean zero used in Theorem 2.1 of [5] to compute the expectations of some 
relevant physical quantities in statistical mechanics. Indeed, their covariances 
coincide since the equality 

cosh (A (t- I)) ^ exp [-At] -b exp [-A (T - t)] 

2Asinh(^) 2A (1 — exp [—AT]) 


holds for every t G [0,T], the period there being the inverse temperature, and 
since the right-hand side of (l87l) is the preferred form to write the covariance in 
[8] and |12j . More generally, certain multidimensional time-periodic processes 
such as T] been useful to solve filtering, smoothing and prediction 

problems as in |14j . and also have played an important role in areas as diverse 
as the analysis of random evolution of loops as in [16) . or the elaboration of 
global shape models generated by periodic lattices as in Chapters 11 and 16 of 
[7], where boundaries of random planar solid objects were investigated. It is our 
contention that the whole family of processes might have an important 

role to play in those areas as well. 

(4) We have already noted that many of the above developments rest on the 
series expansion for Green’s function which we prove in the appendix below. 
We can get expansions of the same kind for Green’s functions associated with 
forward-backward systems of the form 


5tu(x, t) 
u(x, 0) 


and 


-dtv{x,t) 

v{x,T) 


1 / A^ A 

-Axw(x,t)- f Y |x|^-b «:IT(x)j u{x,t), 

iPo(x), 

-- ^Axu(x,t) - ^Y \xf + kW{x)^ v{x,t), 

-- 


(x, t) G X (0, T], 


(x,t) G X [0,T), 


for suitable anharmonic potentials W : i—R where k > 0. We defer 

the analysis of Bernstein processes associated with such systems to a separate 
publication. 
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We complete this article by proving an important property of (1381) which we 
used in the main part of the article. 


5 Appendix: a series expansion for Mehler’s N- 
dimensional kernel 

The notation in this appendix is, of course, the same as in the preceding sections. 
The expansion in question is the following: 


Proposition Al. We have 

9x{x,t,y) 

= ^ (x) X (y) (88) 

mGN" 

where the series converges absolutely for every t G (0,T] and uniformly for all 
X, y S K.^. Furthermore, is indeed Green’s function associated with the 

partial differential equations in m- 


Proof. We first prove the result for the one-dimensional case N = 1, namely. 


f ^ Y 

A (cosh (At) - 1 - — 2xy) 

\27rsinh (At)/ ^ 

2 sinh (At) 


+ 00 

= ^ (a:) (j/) ■ (89) 

m—0 


Our point of departure for this is Mehler’s notable formula in the form 


1 


exp 


TT (1 — 7^) 

+ 00 

= ^ {x) hm {y) 


(l + 7^) (x^ + j/2) - A'yxy 
2 (1-7^) 


771=0 


valid for every 7 G (— 1 ,-|- 1 ), which allows one to express the probability density 
of two jointly Gaussian variables as a power series in the correlation parameter 
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7 (see, e.g., [IS])- By using the scaled Hermite functions (1551) instead, we obtain 


f ^ ^ ^ 

A ((1 + 7 ^) {x"^ + y"^) - ‘^-jxy) 

V^r (1-72)7 

2 ( 1 - 7 ^) 


+ 00 

= X! (x) hm,X iv) , (90) 

m—0 


so that in order to prove (1551) we only need to identify 7. To this end we first 
compare the argument of the exponential on the left-hand side of ((551) with that 
of (1901) . which gives the two conditions 


1-1-7^ 

1-72 


= coth (At) 


1 — 72 2 sinh (At) 

Dividing (15T1) by (1551) then leads to the quadratic equation 

7^ — 2 cosh (At) 7 -I- 1 = 0, 


(91) 

(92) 


whose only solution in the interval (—1, -1-1) is 


7 = 


e 


— Xt 


(93) 


The substitution of (1551) into (150)) then gives (|551) after some simple algebraic 
manipulations. Now, using (1891) we obtain 


e + (x) X (y) 

meN" 

N +CXD 


j—l rrij—O 

A 


2tt sinh (At) 


f N 

]^exp 

i=i 


A (cosh (At) [xj -I- Uj) - 2xjyj) 


2 sinh (At) 


which is (155)) , and the uniform convergence of the series follows from the Cramer- 
Charlier inequality which ensures the uniform boundedness of the h^s in m 
and in their argument (see, e.g.. Section 10.18 in [B] and the references therein). 
Finally, using the series expansion (l88l) and the eigenvalue equation (|56l) of 
Lemma 1 we get 

1 A^ IxP 

5t5A(x,t,y) = -A„5A(x,t,y)-^5A(x,t,y), 

(x,t) e R^x(0,T] 
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for every y S and 


lim gA(x,t,y) = ^ (x) X (y) = ^(x-y) 

meN'^ 

in the sense of distributions since the constitute a complete or¬ 

thonormal system in (R^,C). Consequently, g\ is indeed Green’s function 
associated with the partial differential equation in (IHTl) . ■ 


Remark. Of course, the semi-group composition law 


gx{x,s + t,y) 

/ dzgxi>^,s,z)gxiz,t,y) 
JS.!^ 


(94) 


valid for all x, y G R.^ is inherent in the fact that gx is Green’s function for (ED), 
but follows most directly from the series expansion (l88l) and the orthogonality 
relations 





N 

hrrij (x) X ^ j—ihjij X (x) = • 

J = 1 


Furthermore, we also observe that (1881) brings out the entire pure point spectrum 
of Lemma 1 in the argument of the exponential, a fact that is crucial in our 
construction of Section 4. 
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